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Near-inertial  currents  (NICs)  often  dominate  the  mean  circulation  at  the  East  Flower  Garden  Bank 
(EFGB),  part  of  the  Flower  Garden  Banks  National  Marine  Sanctuary.  The  EFGB,  one  of  several  submerged 
coral  reefs,  is  located  in  the  northwestern  Gulf  of  Mexico,  about  190  km  southeast  of  Galveston,  Texas. 
The  bank  is  about  6  km  wide  in  the  east- west  direction  and  rises  to  within  about  20  m  from  the  surface. 
NICs  near  the  EFGB  are  described  using  current  data  from  5  acoustic  Doppler  current  profilers  that  were 
moored  at  the  edges  of  the  bank  and  on  top  of  the  bank  for  about  a  year.  A  wavelet  analysis  was  used  in 
order  to  better  describe  the  nonstationarity  of  the  NICs.  NICs  were  strongest  during  spring  and  summer 
due  to  their  near  resonant  response  with  sea  breeze  and  the  shallowness  of  the  mixed  layer,  and 
exhibited  a  first-baroclinic-mode  vertical  structure.  NICS  were  generally  larger  near  the  surface  and 
extended  to  the  bottom  on  the  west  side  of  the  EFGB  but  only  to  within  about  20  m  of  the  bottom  on 
the  eastern  side  of  the  bank.  NIC  ellipses  were  nearly  circular  and  rotated  clockwise  above  the  top  of  the 
EFGB  but  became  flatter  and  aligned  with  the  bathymetry  with  increasing  depth;  occasionally,  on  the 
eastern  side  of  the  bank,  the  NIC  vectors  rotated  counterclockwise  due  to  probable  effects  of  lee  vortices 
arising  from  the  mean  flow  interacting  with  the  bank.  Most  energy  input  by  the  wind  at  the  surface  was 
likely  transferred  downward  through  divergence  of  the  meridional  flow  against  the  coastal  boundary. 
The  inertial  currents  were  at  times  more  energetic  than  the  mean  flow,  and  often  accounted  for  more 
than  50%  of  the  total  current  energy. 

Published  by  Elsevier  Ltd. 


1.  Introduction 

Near-inertial  current  (NIC)  oscillations  are  common  to  the  open 
ocean  and  to  shelf  regions  (Webster,  1968).  They  are  ubiquitous  on 
the  shelf  in  the  northern  Gulf  of  Mexico  (Chen  et  al.,  1996; 
DiMarco  et  al.,  2000).  NICs  are  characterized  by  clockwise  (in  the 
northern  hemisphere)  rotations  of  nearly  circular  horizontal  cur¬ 
rents  at  frequencies  at  or  near  the  local  inertial  frequency  (Pollard, 
1970;  Perkins,  1972;  Kundu,  1976).  NICs  are  forced  by  winds  and 
most  effectively  by  wind  pulses  on  time  scales  of  less  than  half  of 
the  local  inertial  period  (Pollard,  1970;  Perkins,  1972).  Existing 
NICs  can  be  enhanced  by  diurnal  sea-breeze  forcing  if  the  currents 
and  sea  breeze  are  in  phase  (near-resonant  conditions  near  30° 
latitude)  (Simpson  et  al.,  2002;  DiMarco  et  al.,  2000;  Jarosz  et  al., 
2007).  Inertial  oscillations  have  been  observed  after  the  passage  of 
hurricanes  (Price,  1976;  Brooks,  1983,  Teague  et  al.,  2007)  and 
atmospheric  cold  fronts  (Chen  et  al,  1996;  Halper  et  al.,  1988; 
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Daddido  et  al.,  1978).  They  can  last  from  several  days  when 
generated  by  local  winds  and  cold  fronts  to  more  than  a  week 
when  generated  by  storms  and  hurricanes.  NICs  can  generate  the 
largest  percentage  of  current  variance,  particularly  where  tides  are 
weak,  for  periods  of  about  a  day.  They  commonly  are  of  the  order 
of  magnitude  of  tidal  variance  in  the  northern  Gulf  of  Mexico 
(Chen  et  al.,  1996). 

Long-term  measurements  of  current  velocity,  temperature,  and 
salinity  were  made  at  the  East  Flower  Garden  Bank  (EFGB)  from 
December  2010  to  December  2011  (Fig.  la).  The  EFGB  is  part  of  a 
marine  sanctuary  designated  as  the  Flower  Garden  Banks  National 
Marine  Sanctuary,  and  is  managed  by  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA).  The  Flower  Gardens,  the 
northern-most  live  coral  reef  in  the  Gulf  of  Mexico,  are  located 
near  the  edge  of  the  Texas-Louisiana  continental  shelf  in  the 
northwestern  Gulf  of  Mexico  (see  http://flowergarden.noaa.gov). 
The  EFGB  is  about  6  km  wide  in  the  along-shelf  direction  and 
10  km  in  length,  and  peaks  at  about  20  m  below  the  surface.  The 
EFGB  is  quite  steep  on  the  eastern  and  southern  sides  and  less  so 
on  the  western  and  northern  sides.  All  of  the  banks  are  sur¬ 
rounded  by  oil  and  gas  activities.  Interactions  of  the  bank  with  the 
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Fig.  1.  (a)  Experiment  area  off  the  Texas-Louisiana  coast  indicated  by  the  box  near 
the  shelf  edge  (100  m  isobath).  Triangles  denote  positions  of  meteorological 
stations,  (b)  Locations  of  ADCP  moorings  (M1-M5)  at  the  East  Flower  Garden  Bank 
(EFGB)  located  near  the  shelf  break  in  the  northwestern  Gulf  of  Mexico.  Bathymetry 
contours  are  every  10  m.  Temperature/salinity  strings  S1-S4  (not  shown)  were 
within  about  200  m  of  corresponding  M1-M4  ADCPs.  The  triangle  denotes  the 
position  of  the  easternmost  meteorological  station. 

currents  are  important  in  processes  over  the  bank.  Protecting  this 
natural  habitat  is  very  important  to  the  ecology  of  the  entire  Gulf 
of  Mexico.  NlCs  are  usually  generated  at  the  surface  by  the  local 
winds,  but  circulation  over  rough  bathymetric  features  at  the  edge 
of  the  shelf,  such  as  the  EFGB,  may  impact  the  structure  and 
dynamics  of  the  NlCs  and  ultimately  influence  submesoscale  mixing 
processes  (Mourn  and  Nash,  2000)  at  this  bank,  which  may  impact 
other  nearby  banks  as  well.  The  NlCs  can  also  contribute  signifi¬ 
cantly  to  the  overall  current  variability  at  the  EFGB. 

The  objective  of  this  work  is  to  characterize  the  impact  of  a 
submerged  bank  near  the  shelf  edge  on  the  current  flow  which  is 
important  for  understanding  of  mixing  and  exchange  mechanisms 
between  the  deep  ocean  and  the  shelf,  and  to  describe  the 
structure  of  near  inertial  currents  over  depth  and  time  at  locations 
surrounding  and  on  top  of  EFGB  under  various  wind  forcing 
conditions.  The  effects  of  banks  on  inertial  oscillations  in  modify¬ 
ing  their  structure  and  dynamics  are  not  now  well  understood  due 


to  a  scarcity  of  current  measurements  next  to  submerged  banks. 
The  structure  of  inertial  oscillations  at  the  EFGB  has  never  been 
studied.  The  EFGB  is  large  enough  to  alter  the  circulation  at  the 
shelf  edge  (Teague  et  al.,  2013)  and  hence  should  have  an  effect  on 
the  inertial  oscillations.  The  NlCs  can  reverse  the  usual  eastward 
mean  flow  over  the  bank  when  their  amplitudes  become  large  but 
these  reversals  may  just  last  a  few  hours  (Teague  et  al.,  2013). 
NlCs  can  impact  the  dispersion  of  contaminants,  biological 
organisms,  and  inorganic  materials  and  hence  may  play  a  key  role 
in  bank  habitats  and  oceanography. 

The  remainder  of  the  paper  is  organized  as  follows.  Current 
data  are  described  in  Section  2,  and  wind  and  pressure  data  are 
described  in  Section  3.  The  methodology  for  calculation  of  the 
mixed-layer  depth  and  a  brief  description  of  the  mixed  layers  are 
provided  in  Section  4.  The  wavelet  analysis  of  the  currents  is 
described  in  Section  5.  Similarly,  a  wavelet  analysis  is  applied  to 
the  winds  in  Section  6.  Results  for  the  diurnal-inertial  band,  based 
on  the  wavelet  analysis,  of  the  currents  are  presented  in  section  7. 
Finally,  some  discussion  and  conclusions  are  given  in  Sections 
8  and  9,  respectively. 


2.  Current  data 

Four  acoustic  Doppler  current  profilers  (ADCPs,  M1-M4)  were 
moored  on  the  bottom  at  the  edge  and  one  ADCP  (M5)  was 
moored  near  the  top  of  the  EFGB  from  December  2010  to 
December  2011.  All  moorings  were  recovered  and  redeployed  in 
June  2011.  Locations  of  the  moorings  and  bathymetry  are  shown  in 
Fig.  lb.  M5  was  located  just  to  the  north  of  the  peak  of  the  bank, 
at  a  depth  of  47  m.  M1-M4  were  deployed  just  NW,  SW,  SE  and 
NE,  respectively,  of  the  EFGB,  at  depths  ranging  from  100  to  127  m. 

RD  Instruments  Workhorse  ADCPs,  operating  at  300  kH  (Ml- 
M4)  and  600  kHz  (M5),  were  mounted  in  trawl-resistant  pods  that 
rested  on  the  bottom  and  are  referred  to  as  Barnys  due  to  their 
barnacle-like  shape  (Perkins  et  al.,  2000).  They  recorded  nearly  full 
water  column  current  profiles  at  2  m  vertical  resolution  except  for 
4  m  vertical  resolution  at  M3  every  12  or  15  min.  ADCP  accuracy  is 
0.5%  of  the  water  velocity.  The  usual  naming  conventions  for  the 
velocities  are  used:  U  positive  towards  the  east,  V  positive  towards 
the  north,  and  W  positive  upwards.  Further  details  on  the  data  can 
be  found  in  Teague  et  al.  (2013).  After  initial  processing  and 
editing,  ADCP  data  were  resampled  at  even  30-min  times  (com¬ 
patible  with  wind  data),  and  time-trimmed  (at  the  series  ends)  so 
that  all  series  from  deployment  1  (Dl)  started  and  stopped  at  the 
same  times  (except  for  M5  which  was  short),  and  similarly  for 
deployment  2  (D2). 

Above  the  EFGB  (shallower  than  about  40  m)  horizontal  velo¬ 
cities  are  significantly  correlated.  EOF  analyses  applied  to  U  and  V 
at  20  m  and  40  m  depths  show  that  more  than  80%  of  overall 
variance  was  in  the  first  horizontal  EOF  mode  and  less  than  10%  in 
mode  2  at  20  m.  The  mode  1  amplitudes  are  virtually  constant  for 
all  moorings  used  (M1-M4  in  deployment  1,  M1-M5  in  deploy¬ 
ment  2),  and  mode  1  explained  variance  for  each  location  is  at 
least  75%  at  20  m.  Overall  first-mode  variance  decreases  with 
depth  with  an  explained  variance  minimum  of  60%  at  40  m.  This 
result  implies  that  on  the  relatively-small  horizontal  extent  of  the 
observations  that  currents  above  the  depth  of  the  EFGB  are 
laterally  highly  coherent  and  in  phase. 

Inertial  periods  at  the  center,  northern,  and  southern  edges  of 
the  EFGB,  using  the  100  m  isobath  as  the  boundary,  are  approxi¬ 
mately  25.55  h,  25.52  h,  and  25.58  h,  respectively.  The  3  major 
tidal  components  are  M2,  Kl,  and  01  and  have  periods  of  12.42  h, 
23.93  h,  and  25.82  h,  respectively.  The  ADCP  data  were  harmoni¬ 
cally  analyzed  for  4  diurnal  and  4  semidiurnal  tidal  constituents 
(l<i,  Oi,  Pi,  Qi,  N2,  M2,  S2, 1<2).  This  analysis  was  done  using  depth- 
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averaged  currents  at  1VI1-M4;  D1  and  D2  records  were  concate¬ 
nated  yielding  overall  record  lengths  of  nearly  one  year  (363.8 
days).  The  MS  records  were  excluded  since  1VI5-D1  was  only  40 
days  long.  The  Rayleigh  criteria  for  most  combinations  of  the 
selected  constituents  are  about  27  days  or  less;  for  Ki  and  Pi,  as 
well  as  for  1<2  and  S2,  it  is  182.6,  well  below  the  record  length  of 
363  days;  however,  the  year-long  record  is  necessary  for  1<2  due  to 
possible  aliasing  with  a  semiannual  constituent.  There  were  small 
variations  of  amplitudes  among  the  mooring  sites  which  may  be 
due  to  shoaling  effects  near  the  shelf  edge  and  the  bank.  The 
amplitudes  for  Oi,  l<i  and  M2  were  approximately  1.9cms“\ 
2.3cms“^  and  2.2cms“\  respectively.  The  other  5  constituent 
amplitudes  were  small  (  <  1  cm  s“^).  These  amplitudes  are  some¬ 
what  smaller  than  those  reported  by  DiMarco  and  Reid,  (1998). 
Amplitudes  of  the  inertial  oscillations  can  exceed  15cms“^  and 
are  often  larger  than  tidal  amplitudes.  NlCs  can  be  modified  by  the 
diurnal  tide,  depending  on  the  phase  difference,  but  this  effect  is 
relatively  small.  Rotary  autospectra  for  current  velocity  at  7.5  m 
depth  from  Ml  are  displayed  in  Fig.  2  and  show  a  strong  clockwise 
(CW)  peak  near  the  diurnal-inertial  as  well  as  semidiurnal 
periods.  There  is  little  energy  at  the  diurnal  frequency  for  the 
counterclockwise  (CCW)  component  but  some  CCW  semi-diurnal 


Fig.  2.  Rotary  autospectra  for  ADCP  M3-D2.  CW:  with  tides  (black);  barotropic 
tides  removed  (green).  CCW:  with  tides  (red);  barotropic  tides  removed  (blue). 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


energy  is  present  but  at  a  lower  level  than  for  the  CW  component. 
Removal  of  principal  barotropic  tidal  constituents  (also  Fig.  2) 
attenuates  much  of  the  semidiurnal  energy  but  the  CW  diurnal- 
inertial  peak  is  minimally  affected. 

Overall  average  squared  vertical  shear  of  horizontal  currents 
(5h^  =  {dUldzf -TidVIdz)^)  ranged  from  0.005  s“^  to  0.011  s“^  at 
mooring  sites.  The  squared  shear  combined  with  vertical  density 
gradients  (based  on  temperature  and  salinity  data  from  string 
moorings.  Section  4)  were  used  to  calculate  gradient  Richardson 
numbers,  R/  =  N^/Sh^,  where  is  the  buoyancy  frequency.  With  a 
few  exceptions,  some  of  which  were  likely  spurious,  R/  values  were 
above  0.25  indicating  vertical  stability  over  most  of  the  record. 
Ri  values  in  the  mixed  layer  were  not  included  since  both  vertical 
shear  and  density  gradients  were  small. 


3.  Wind  and  pressure  data 

Wind  velocity  and  atmospheric  pressure  data  were  collected 
half-hourly  at  the  southern  edge  of  the  EFGB  at  National  Data  Buoy 
Center  Station  42047,  and  at  Station  42046  located  approximately 
37  km  west  of  the  EFGB.  Winds  were  very  similar  between  these 
two  stations  for  the  concurrent  recording  periods.  Winds  at  42046 
were  used  when  42047  was  out  of  service  which  was  about  half  of 
the  year-long  mooring  period.  Winds  during  the  deployment 
period  were  typical  for  this  region  (de  Velasco  and  Winant, 
1996).  Wind  velocities  and  pressures  are  shown  in  Fig.  3.  Wind 
gusts  are  considerably  larger  than  the  30-min  averaged  velocities. 
The  atmospheric  pressure  record  shows  numerous  frontal  pas¬ 
sages  from  December  2010  to  May  (days  -30  to  130)  and  from 
October  2011  through  December  2011  (days  270-360).  Winds 
were  weakest  mainly  during  summer  (except  for  a  wind  event 
near  day  245),  and  were  generally  southerly.  Hurricanes  or  tropical 
storms  did  not  impact  the  EFGB  in  2011  and  atmospheric  pressure 
was  almost  always  above  1000  mb  (Fig.  3).  There  is  a  significant 
peak  in  the  clockwise  (CW)  rotary  wind  autospectrum  (Fig.  4)  for 
the  summer  period  that  may  be  associated  with  sea-breeze  effects 
(Zhang  et  al.,  2009).  The  most  intense  inertial  oscillations  were 
found  during  spring  and  summer.  These  enhanced  amplitudes  are 
associated  with  a  near-resonant  response  to  sea-breeze  effects 
(Teague  et  al.,  2013;  Zhang  et  al.,  2009;  Jarosz  et  al.,  2007;  Simpson 
et  al.,  2002;  DiMarco  et  al.,  2000).  The  EFGB  is  about  200  km  off 
the  Texas-Louisiana  coastline  and  is  considerably  farther  offshore 


Time  (days  in  2011) 


Fig.  3.  Eastward  (a)  and  northward  (b)  components  of  wind  velocity;  atmospheric  pressure  (c),  near  the  EFGB.  First  sample  on  1  December  2010. 
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Frequency  (cpd) 


Fig.  4.  Clockwise  (CW)  rotary  autospectra  of  the  wind  velocity  for  each  season 
(black-winter,  red-spring,  summer-green,  fall-blue)  are  shown.  A  significant  diurnal 
peak  is  present  for  the  summer  (July-September).  CCW  spectra  (not  shown)  are 
similar,  but  a  diurnal  peak  is  not  present.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 

than  the  50  km  range  that  Chen  et  al.  (1996)  suggested  that  sea- 
breeze  effects  could  be  significant.  However,  Zhang  et  al.  (2009) 
found  that  the  sea  breeze  can  extend  to  300  km  off  the  coast  in 
this  region. 


4.  Mixed-layer  depths 

Temperature  data  from  string  moorings  located  near  (within 
several  hundred  meters)  M1-1VI4  were  used  to  estimate  the 
mixed-layer  depth  (MLD)  (Teague  et  al.,  2013).  Temperature  from 
string  mooring  instruments  (typically  10)  were  interpolated  to 
one-meter  intervals,  but  not  extrapolated  to  the  surface.  The  MLD 
was  defined  as  the  depth  where  the  temperature  changed  by  at 
least  0.25  °C  from  the  shallowest  temperature  measurement.  The 
MLDs  computed  from  each  string  were  smoothed  over  24  h  to 
reduce  the  higher  frequency  variability.  Finally,  the  average  MLD 
over  all  the  strings  was  computed  and  a  five-day  boxcar  average 
was  applied  to  reduce  the  daily  variability.  MLDs  were  deepest, 
60-80  m,  during  December  and  January.  Throughout  much  of  the 
year,  MLDS  ranged  between  20  and  40  m.  Rapid  changes  in  the 
MLD  may  have  been  caused  by  eddy  passages,  such  as  the  rapid 
decrease  in  MLD  during  March  (Teague  et  al.,  2013).  The  seasonal 
cycle  in  MLD  at  the  EFGB  is  similar  to  that  reported  by  McGrail 
(1983).  The  average  MLD  is  shown  in  the  upcoming  wavelet  plots 
(Figs.  12-19). 


5.  Wavelet  analysis 

Currents  at  the  EFGB  change  quickly  over  time  due  to  passing 
eddies,  rapidly  changing  wind  fields,  and  interaction  with  the  bank. 
Hence  the  current  patterns  are  not  very  predictable  (Teague  et  al., 
2013).  These  currents  are  not  stationary  and  conventional  Fourier 
spectrum  analysis,  which  provides  information  about  the  average 
amplitude  and  phase  for  each  harmonic  and  energy  content  for  the 
entire  time  span  of  the  data,  presumes  the  data  are  stationary.  This 
method,  however,  does  not  provide  information  about  the  energy 
levels  and  its  variations  within  the  time  span  of  the  data  records. 
Wavelet  analyses  of  ocean  currents  (Liu  and  Miller,  1996)  can 
provide  an  effective  way  of  obtaining  time-frequency  information 
of  non-stationary  processes  such  as  NlCs.  Wavelet  analysis  provides 
a  complementary  approach  to  the  traditional  Fourier  spectrum 
analysis,  and  is  used  here  to  provide  a  description  of  the  time 
evolution  of  NlCs  near  the  EFGB.  Significant  current  fluctuation 
energy  is  present  in  the  ADCP  velocity  series  acquired  near  the 


EFGB.  The  objectives  here  are  to  show  the  importance  of  variability 
in  the  diurnal-inertial  band  (DIB)  relative  to  fluctuations  at  other 
frequencies,  to  describe  the  characteristics  of  variability  in  the  DIB, 
to  examine  the  relationship  of  DIB  variability  with  local  hydro- 
graphic  conditions,  and  to  calculate  coherences  between  currents 
and  wind  velocity  over  the  area. 

Wavelet  decomposition  of  a  time  series  is  essentially  the 
systematic  application  of  a  set  of  band-pass  Alters  whose  band- 
widths  are  proportional  to  the  center  periods  of  each  (Torrence 
and  Compo,  1998;  Mallat,  1989).  Here,  we  apply  continuous 
wavelet  transforms  (CWTs)  using  a  complex  Morlet  wavelet,  with 
MATLAB  software  developed  by  Torrence  and  Compo  (1998). 
The  complex  Morlet  wavelet  is  essentially  a  cosine,  sine  wave 
with  a  Gaussian  envelope  for  the  real,  imaginary  parts,  respec¬ 
tively;  these  functions  are  then  convolved  with  the  data  to 
produce  a  set  of  band-passed  complex  series  for  each  variable 
where  the  imaginary  part  is  shifted  by  -90°  from  the  real  part. 
Wavelet  decomposition  was  done  with  the  half-hourly-sampled 
current  series,  and  with  the  local  wind  data. 

The  wavelet  analysis  results  in  complex  series  of  U  and  V 
velocity  components  which  are  functions  of  depth,  time  and 
frequency  band,  e.g. 

u(z,  t)  =  i;  u(z,  tjfc)  =  i;  (Ur(z,  tjfc) + iUiiz,  tji,)) ,  (i  > 

where  the  tilde  indicates  wavelet  decomposition  and  k  corre¬ 
sponds  to  a  particular  band.  Following  Gonella  (1972),  and  Mooers 
(1973)  band-limited  velocities  can  be  transformed  to  rotary 
components 

20^  =0r-Vi, 

2V+  =Ui  +  Vr, 

2U~  =  0r  +  Vi, 

2V~  =Ui-Vr,  (2) 

where  the  plus  and  minus  superscripts  indicate  CCW  and  CW 
rotation,  respectively,  of  the  velocity  vector.  We  calculate  specific 
energy 

S=|&|V|v|^  (3) 

for  total  velocity;  energies  for  CCW  and  CW  currents,  5^  and  S  , 
are  similarly  defined.  Plots  of  S  and  for  depths  of  20  m  (Ml- 
M5),  40  m  (for  M5)  and  70  m  (M1-M4)  appear  as  Figs.  5-8. 
The  y-axis  shows  the  period  instead  of  the  frequency  for  easier 
visualization.  The  inertial  oscillations  will  appear  in  the  CW 
spectra  at  periods  just  greater  than  a  day  for  the  latitude  of  the 
EFGB.  Because  the  inertial  oscillations  rotate  clockwise,  they  are 
not  evident  in  the  CCW  spectra  in  the  Northern  Hemisphere. 

The  most  prominent  features  in  Fig.  5  ( ~20  m)  are  the  bursts  of 
CW  DIB  energy,  strongest  between  days  100  and  250,  at  all 
mooring  locations.  This  should  not  be  surprising  given  the  diurnal 
peak,  associated  with  the  sea  breeze,  in  the  CW  wind  spectrum 
(Fig.  4)  during  this  time  period.  The  high  CW  energy  near  the 
inertial  period  at  day  245  is  associated  with  a  large  north- 
northwesterly  wind  burst  (Fig.  3).  Some  bursts  of  semi-diurnal 
energy  are  also  evident,  especially  near  day  170  but,  as  indicated 
above,  much  of  this  may  be  tidal.  Lower-frequency  (LF)  energy  is 
also  present,  but  is  generally  lower  than  that  in  the  DIB.  In 
particular,  there  are  several  bursts  at  periods  of  several  days  near 
times  250,  260  and  320  days.  Also,  between  days  220  and  280 
there  is  a  persistent  burst  of  ~  17-day-period  energy,  primarily  at 
M2,  M3  and  M5.  DIB  energy  is  essentially  absent  in  plots  of  CCW 
energy  (Fig.  6).  The  LF  features,  on  the  other  hand,  are  more 
pronounced  and  also  show  up  prominently  at  all  mooring  loca¬ 
tions,  with  quite  similar  signatures,  between  days  230  and  280. 
LF  CCW  features  also  appear  near  days  40-50,  primarily  at  M2  and 
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Fig.  5.  Clockwise  (CW)  energy  for  currents  near  depth =20  m,  as  functions  of  period  and  time.  Green  areas  indicate  truncation  limits,  most  evident  for  longer  periods. 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  6.  Counterclockwise  (CCW)  energy  for  currents  near  depth =20  m,  as  functions  of  period  and  time.  Green  areas  indicate  truncation  limits,  most  evident  for  longer 
periods.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  7.  Clockwise  (CW)  energy  for  currents  near  depths  70  m  (M1-M4)  and  40  m  (M5),  as  functions  of  period  and  time.  Green  areas  indicate  truncation  limits,  most  evident 
for  longer  periods.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  8.  Counterclockwise  (CCW)  energy  for  currents  near  depth=70  m  (M1-M4)  and  40  m  (M5),  as  functions  of  period  and  time.  Green  areas  indicate  truncation  limits,  most 
evident  for  longer  periods.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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M3.  These  LF  motions  are  likely  associated  with  passing  eddies,  as 
suggested  during  these  time  periods  by  Teague  et  al.  (2013).  CW 
energies  near  70  m  (M1-M4)  and  40  m  (M5)  (Fig.  7)  are  quite 
similar  to  the  shallower  case  (Fig.  5),  but  show  some  attenuation  of 
the  DIB  levels.  Much  of  the  ~  17-day  energy  apparent  at  20  m  is 
gone  at  70  m,  but  at  Ml  there  is  a  significant  burst  of  several-day 
energy  near  day  245  which  also  corresponds  to  the  wind  burst  at 
the  end  of  summer.  Deeper  CCW  energies  (Fig.  8)  also  show 
reduction  in  DIB  intensity,  but  not  as  much  as  for  the  shallower 
level.  Note  the  significant  deep  CCW  DIB  energy  bursts  at  M2  and 
particularly  M3  near  day  245;  these  may  be  related  to  bathymetric 
effects  of  the  EFGB. 


6.  Wind  forcing 

Wavelet  analysis  also  was  applied  to  the  wind  velocities. 
Period-time  plots  of  CW  and  CCW  (Sw)  appear  in  Fig.  9.  Seasonal 
changes  in  the  energy  are  evident  with  broadband  frontal  pas¬ 
sages,  primarily  from  late  fall  to  early  spring.  Emergence  and 
intensification  of  diurnal  energy,  mostly  CW,  are  found  during  late 
spring  and  summer  into  late  fall.  Between  days  170  and  270,  the 
peak  in  energy  for  the  DIB  is  well  defined  and  is  indicative  of  the 
strong  sea  breeze.  Strong  broadband  summer  events  occur  near 


days  170  and  210.  There  is  some  CW  energy  in  the  semidiurnal 
band  during  summer;  it  may  be  present  at  other  times  but  is 
obscured  by  frontal  passages  and  storms.  We  investigate  possible 
wind  forcing  of  the  observed  NlCs  with  a  wavelet  coherence 
algorithm  (Grinsted  et  al.,  2004;  Torrence  and  Compo,  1998). 
The  cross-wavelet  transform  (Cww)  for  northward  wind  and 
current  (at  ~20  m)  is  defined  by 

(4) 


where  the  asterisk  denotes  complex  conjugate.  To  properly  calcu¬ 
late  coherence  {yww),  smoothing  over  time  and  frequency 
(Grinsted  et  al.,  2004),  indicated  by  an  overbar,  is  applied  to  C, 

I  ~  I  ~  |2 

V  ,  and  Vw  .  and 


7^  vvw(f’  0  — 


\CvVw(f ,  t)\ 


(5) 


The  phase  of  the  current  relative  to  wind  fluctuations  is  given 
by 

^vvwd,  0  =  tan  -^CvvwCf,  t),  (6) 


Coherences  of  noise  series  (using  a  Monte  Carlo  Method)  yield 
a  5%  significance  level,  Y5(f)  (only  5%  of  coherences  with  random 
data  rise  above  this  level).  This  provides  a  measure  of  confidence 
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Fig.  9.  CW  and  CCW  wavelet  wind  energies.  The  time  series  was  split  to  match  the  time  intervals  of  the  two  ADCP  deployments.  Green  areas  indicate  truncation  limits,  most 
evident  for  longer  periods.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  10.  Coherence-squared  between  ~20  m  northward  current  and  northward  wind  at  Ml.  The  5%  significance  level  is  about  0.85  for  periods  between  2  h  and  about  50 
days,  increasing  gradually  for  longer  periods.  Green  areas  indicate  truncation  limits,  most  evident  for  longer  periods.  24-h  period  is  indicated  by  the  black  line.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  11.  Temporal  phase  of  ~20  m  northward  current  relative  to  northward  wind  at  Ml,  for  3  adjacent  diurnal-inertial-band  periods  19.656  h,  23.375  h,  27.798  h;  positive 
phase  indicates  current  leading  wind. 
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for  Yww  75  is  a  function  of  frequency  and  ranges  from  0.84  to  0.86 
between  periods  of  1.5  h  and  40  days,  For  longer  periods  it 
increases  nearly  linearly,  exceeding  0.90  at  60  days,  ys  rises  sharply 
for  periods  less  than  1.5  h. 

Fig.  10  shows  /vvw  at  Ml ;  the  other  locations  are  qualitatively 
similar,  as  are  plots  of  yuuw  (not  shown).  The  northward  compo¬ 
nent  was  selected  here  to  align  with  the  principal  direction 
(onshore-offshore)  of  the  diurnal  sea  breeze,  which  is  an  impor¬ 
tant  forcing  agent.  The  similarities  among  locations  should  not  be 
surprising  since  correlations  among  20  m  current  velocities  from 
the  moorings  are  quite  high.  Significant  coherences  (orange  to  red 
to  brown)  are  most  evident  in  the  diurnal-inertial  range  and  at 
several-day  to  several-week  periods.  The  latter  correspond  to 
energetic  meridional  wind  events,  e.g.  near  day  250  (see  Figs.  9, 
3b).  Phases  {0ww  )  for  the  3  periods  included  in  the  DIB  (Fig.  11), 
plotted  only  where  yww  is  significant,  show  periods  of  reinforcing 
or  resonance  ( -90°  <  0  <  90°),  particularly  between  days  120  and 
270.  There  are  also  a  few  short  periods  consistent  with  attenuation 
(101  >  90°),  e.g.  near  day  230. 


7.  Diurnal-inertial  band 

Zhang  et  al.  (2009)  define  the  band  between  periods  of  20  h 
and  28  h  as  the  diurnal  inertial  band  (DIB).  In  this  analysis,  the 
wavelet  band  central  periods  (in  hours)  include  the  following:  [..., 
16.5290, 19.6560,  23.3750,  27.7980,  33.0570,. . .  ]  (out  of  a  total  of  53 
bands),  so  we  define  our  DIB  as  the  2nd-4th  of  the  5  periods 
shown.  Flence,  we  sum  currents  for  these  three  bands  to  form  DIB 
velocity 

20  ^ 

UD(Z,t)=  Z  U(z,t,f0;  (7) 

/<  =  18 


likewise  for  Vd  and  corresponding  rotary  components.  It  is 
useful  to  examine  some  time-dependent  statistics  of  the  DIB 
velocities.  In  the  following,  the  statistics  are  50%-overlapping 
two-day  boxcar  averages  of  respective  quantities,  yielding  daily 
values.  The  time-smoothed  quantities  are  denoted  by  overbars. 

The  smoothed  energies.  So,  ,  can  be  interpreted  as  time 

and  depth  dependent  spectra  for  the  DIB.  The  CW  spectral 
series,S^  (z,  t),  is  shown  in  Fig.  12.  The  signature  of  5^  is  nearly 
the  same  as  that  of  So  (not  shown),  indicating  that  CW  energy 
dominates.  The  CCW  spectra,  5^  (z,  t)  (Fig.  13),  show  a  much  lower 
energy  throughout  most  of  the  record.  Notable  exceptions  include 
bursts  near  70  m  and  day  170  at  M2  and  M3,  and  between  day  140 
and  230  at  M4.  There  are  other  CCW  bursts  of  shorter  extent  in 
time  and  depth.  Where  the  CCW  and  CW  energy  are  comparable, 
the  currents  are  approximately  rectilinear.  The  highest  CW 
energy  levels  (Fig.  12)  occur  primarily  during  the  spring  and 
summer  (days  110-240),  but  there  are  other,  more  isolated,  bursts 
outside  this. 

In  general,  the  current  vectors  of  band-limited  rotary  motions 
describe  ellipses.  These  ellipses  can  range  from  circular  to  recti¬ 
linear  in  shape  with  either  CCW  or  CW  rotation.  Parameters  of  the 
ellipses,  which  here  are  functions  of  depth  and  time,  can  be 
derived  from  the  spectra  and  other  statistical  series.  The  semi¬ 
major  and  semi-minor  axes  are  given  by 


A  positive  Xmm  indicates  CCW  rotation;  Xmin  =  0  results  if  the 
motion  is  rectilinear.  If  Xmin=  ±^maj  the  motion  is  circular.  When 
the  magnitude  of  the  ratio  R=Xminl^maj  is  significantly  less  than 
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Fig.  12.  Clockwise  (CW)  spectra  for  diurnal-inertial-band  (DIB)  currents  as  functions  of  depth  and  time.  Statistics  averaged  over  2-day  segments.  The  black  curve  represents 
mixed-layer  depth. 
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Fig.  13.  Counterclockwise  (CCW)  spectra  for  diurnal-inertial-band  (DIB)  currents  as  functions  of  depth  and  time.  Statistics  averaged  over  2-day  segments.  The  black  curve 
represents  mixed-layer  depth. 
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Fig.  14.  Semi-major  axis  length  for  DIB  current  ellipses.  Green  areas  indicate  values  in  excess  of  25  cm  s  \  The  black  curve  represents  mixed-layer  depth.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  15.  Ratio  of  minor  to  major  axis  for  DIB  current  ellipses.  Negative  values  indicate  clockwise  rotation  of  the  current  vector.  The  black  curve  represents  mixed-layer  depth. 
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Fig.  16.  Orientations  of  DIB  current  ellipses.  Positive  values  are  counterclockwise  from  east-west.  Green  areas  indicate  orientations  that  are  indeterminate  due  to 
uncertainties.  The  black  curve  represents  mixed-layer  depth.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
this  article.) 
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Fig.  17.  Ratio  of  DIB  energy  to  total  energy.  The  green  curve  represents  mixed-layer  depth.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 
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Fig.  18.  DIB  zero-lag  correlation  for  V2)(z,t)  relative  to  VoiZoT),  where  Zq  is  the  depth  of  the  shallowest  bin.  The  green  curve  represents  mixed-layer  depth. 
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unity,  it  is  useful  to  calculate  ellipse  orientation 


xif  =  tan  ^ 


2UdVd  \ 


(9) 


Uncertainty  in  y/  increases  as  R  approaches  +  1.  X^q/  (plotted  in 
Fig.  14)  represents  the  smoothed  DIB  current  amplitude;  the 
information  here  is  similar  to  that  shown  in  Fig.  12,  but  is  more 
directly  comparable  to  current  velocity.  The  ellipse  axis  ratios  (R) 
(Fig.  15)  show  the  dominance  of  CW  motions  in  the  DIB  and 
emphasize  the  CCW  periods,  especially  at  M4  deeper  than  the  top 
of  the  EFGB.  The  occurrences  of  CCW  energy  at  M4,  which  was 
located  near  the  steep  northeast  side  of  the  EFGB,  may  be 
associated  with  the  generation  of  cyclonic  vortices  in  the  lee  of 
the  bank  when  mean  currents  are  eastward  Qarosz  et  al.,  submitted 
for  publication).  These  CCW  periods  in  DIB  energy  roughly  coincide 
with  periods  of  eastward  mean  flow  reported  by  Teague  et  al.  (2013; 
their  Fig.  8a).  At  all  locations,  lighter  colors  (smaller  magnitude  of 
the  ratio)  indicate  flatter  ellipses,  or  more  rectilinear  fluctuations. 
These  features  suggest  bathymetric  effects  on  the  flow.  Ellipse 
orientations  (or  principal  directions  of  variability)  appear  in 
Fig.  16.  Much  of  the  plot  shows  green  areas;  these  are  where  the 
ratio  \R\  ~  1  and  the  uncertainty  in  calculated  y/  is  greater  than  about 
30°,  rendering  it  meaningless.  At  locations  where  useable  orienta¬ 
tions  exist  (deeper  than  about  60  m),  they  are  about  30°-45°  at  Ml 
and  M3  and  -45°  to  -80°  at  M2  and  M4;  these  values  align 
roughly  with  the  bathymetry  of  the  EFGB  at  these  depths. 

We  can  examine  the  ratio  of  energy  in  the  DIB  relative  to  total 
energy. 


Re{z,  t)  = 


(10) 


where  the  overbar,  as  above,  represents  a  2-day,  overlapping  time 
average.  Re  plotted  in  Fig.  17  is  generally  highest  during  spring- 
summer,  with  the  most  pronounced  event  occurring  between 
days  160  and  180.  The  DIB  energy  often  exceeded  50%  of  the  total 
current  energy  and  sometimes  even  exceeded  70%  of  the  total 
energy.  These  high  energy  contributions  often  extended  from  near 
surface  to  near  bottom. 

Complex  correlations  (recall  that  Ud  and  Vq  are  complex  series) 
calculated  (over  each  2-day  segment,  as  above)  between  velocities 
at  each  depth  relative  to  the  shallowest  depth  yield  series  of 
magnitudes  {p)  and  phases  (0).  The  product  p  cos{0)  yields  zero- 
lag  values  for  Ud  and  Vq.  The  zero-lag  correlation  for  is  shown 
in  Fig.  18;  the  zero-lag  correlation  for  Ud  is  quite  similar  with  the 
main  differences  occurring  where  U  values  are  lower,  particularly 
at  M4  below  the  depth  of  the  top  of  the  EFGB.  The  correlations  are 
concentrated  near  1  or  -1,  suggesting  first-mode  baroclinic 
motion.  Also,  estimates  of  mixed-layer  depths  (indicated  on 
Fig.  18)  correspond  roughly  to  the  transition  from  positive  to 
negative  correlations.  Phases  (0)  are  shown  in  Fig.  19.  Since  these 
results  are  quite  noisy,  the  phases  are  assigned  quadrants  to 
provide  estimates  of  gross  changes  over  depth  and  time.  In  phase 
(315°-45°),  leading  (45°-135°),  opposite  phase  (135°-225°)  and 
lagging  (225°-315°)  are  represented  by  yellow,  red,  green,  blue, 
respectively.  Above  the  mixed-layer  depth  the  current  is  roughly  in 
phase  with  that  at  the  shallowest  depth.  Below  this,  the  current  is 
primarily  oppositely  directed.  There  are  periods  of  leading  and 
lagging  fluctuations  suggesting  upward  and  downward  phase 
propagation  (or  downward  and  upward  energy  propagation), 
respectively,  but  overall,  the  vertical  structure  of  DIB  variability 
suggests  a  first  baroclinic  mode.  Phases  and  hence  zero-lag 
correlations  are  questionable  where  DIB  amplitudes  are  low.  The 
deep  Angers  near  (e.g.)  days  235  and  265  at  M4-D2  (Fig.  18) 
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Fig.  19.  Lagged  correlation  phase  (0)  for  VoizT)  relative  to  VuiZoT),  where  Zq  is  the  depth  of  the  shallowest  bin.  Phases  are  expressed  in  quadrants  centered  on  0°,  90°,  180° 
and  270°.  The  black  curve  represents  mixed-layer  depth.  Currents  are  in  roughly  in  phase  for  the  -45°-45°  interval  (yellow),  have  upward  propagating  phase  for  the  45°- 
135°  interval  (red),  have  opposite  phase  for  the  135°-225°  interval  (green),  and  have  downward  propagating  phase  for  the  225°-315°  interval  (blue).  (For  interpretation  of 
the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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correspond  to  velocity  amplitudes  of  about  5  cm  s  ^  or  less  and 
are  primarily  due  to  uncertainties  in  phase. 


8.  Discussion 

Near-inertial  currents  are  an  important  energy  exchange 
mechanism  between  the  winds  and  ocean  and  can  be  responsible 
for  a  substantial  fraction  of  the  kinetic  energy  throughout  the 
water  column.  Our  findings  in  the  preceding  sections,  as  well  as 
the  high  horizontal  coherence  of  currents  above  the  bathymetry 
(Section  2),  imply  that  NlCs  in  the  vicinity  of  the  EFGB  are 
primarily  forced  by  winds  acting  on  the  mixed  layer.  Fluctuations 
at  or  near  the  inertial  period  (here,  about  25.5  h)  result  from 
broadband  excitation  due  to  storms  and  frontal  passages  (Chen 
and  Xie,  1997),  while  near-resonant  diurnal  (24  h)  responses  often 
occur  when  a  significant  sea  breeze  is  present,  mainly  during  late 
spring  and  summer  months  (DiMarco  et  al.,  2000;  Jarosz  et  al., 
2007;  Simpson  et  al,  2002,  Rippeth  et  al.,  2001).  Furthermore,  the 
response  to  wind  forcing  is  enhanced  when  the  mixed  layer  is 
shallow  and  stratification  is  stronger.  DIB  energy  is  input  to  the 
mixed  layer  as  primarily  CW  currents  with  nearly  circular  ellipses. 
There  is  some  evidence  here  of  energy  generated  near  the  surface 
being  propagated  downward  and  energy  generated  near  the 
bottom  (by  bathymetric  interaction)  being  propagated  upward 
(Kundu,  1976;  Chant,  2001;  Lerczak  et  al.,  2001).  However,  the 
dominant  mechanism  in  this  case  is  likely  due  to  the  presence  of 
the  east-west  boundary  at  the  coast.  This  boundary  allows  a 
meridional  pressure  gradient  to  develop,  associated  with  the 
north-south  component  of  the  NIC,  which  results  in  an 
oppositely-directed  current  below  the  pycnocline  yielding  a  first 
baroclinic  modal  structure  (Rippeth  et  al.,  2001;  Millot  and 
Crepon,  1981;  Pettigrew,  1981;  Kundu  et  al.,  1983;  Tintore  et  al., 
1995;  Chen  et  al.,  1996). 

DIB  energy  penetrates  most  of  the  water  column  rather  quickly 
(within  about  1  day).  Above  the  top  of  the  EFGB  the  NIC  current 
vectors  rotate  CW  and  their  ellipses  are  nearly  circular.  Below  the 
mixed  layer  the  velocities  become  anticorrelated  with  the  mixed- 
layer  velocities.  The  ellipses  flatten  as  depths  increase  beyond 
about  50  m  with  their  major  axes  aligning  roughly  with  the  nearby 
bathymetry,  implying  significant  interaction  with  the  bank.  The 
interaction  can  cause  the  flow  to  become  rectilinear.  In  some  cases 
(mostly  at  1VI4)  the  ellipses  exhibit  the  opposite  vector  rotation 
(CCW)  for  extended  periods;  this  occurs  primarily  when  mean 
flow  (see  Fig.  8a  in  Teague  et  al.,  2013)  is  eastward  over  much  of 
the  water  column  at  Ml  and  M2,  during  spring  and  summer,  when 
DIB  energy  is  highest.  The  lee  side  of  the  EFGB  is  very  steep  while 
the  west  side  slope  is  much  less  extreme.  The  eastward  flow  does 
not  extend  as  deep  at  M3  and  M4  due  to  the  blocking  effect  of  the 
EFGB.  However,  the  presence  of  the  bank  in  the  eastward  flow  can 
produce  cyclonic  vortices  (Jarosz  et  al.,  submitted  for  publication) 
on  the  lee  side  of  the  bank  which  may  further  interact  with  the 
NlCs  resulting  in  the  CCW  rotations  seen  at  M4. 


9.  Conclusions 

NlCs  at  the  EFGB  occur  in  bursts  throughout  the  year  but  are 
most  concentrated  during  the  spring  and  summer  periods  where 
they  are  enhanced  by  near-resonant  sea-breeze  effects.  The  bank 
tends  to  modify  NIC  characteristics  below  the  top  of  the  bank.  NlCs 
can  occasionally  reverse  the  predominantly  eastward  mean  flow 
on  short  time  scales.  Typically,  the  vertical  structure  of  the  NlCs 
followed  a  flrst-baroclinic  mode  with  a  180°  phase  shift  between 
upper  and  lower  layers  but  was  more  complicated  when  waters 
were  weakly  stratified.  NlCs  were  generally  stronger  at  the  surface 


but  sometimes  were  largest  at  mid-depth.  Barotropic  tidal  cur¬ 
rents  were  small  and  generally  less  than  the  NlCs.  Above  the  bank, 
NlCs  were  highly  coherent.  Below  the  top  of  the  bank,  the  NIC 
ellipses  tended  to  align  with  the  bathymetry  and  become  more 
rectilinear  while  flowing  around  the  bank.  CCW  rotation  of  the 
ellipse  vector  was  common  on  the  lee  or  eastern  side  of  the  bank 
below  the  top.  NlCs  were  often  a  significant  part  of  the  total 
current  and  could  account  for  more  than  50%  of  the  total  current 
energy  throughout  the  water  column.  NlCs  are  likely  to  be  partially 
responsible  for  enhanced  mixing  of  physical,  biological,  and 
geological  properties. 

The  Flower  Garden  Banks  are  very  important  to  the  ecosystem 
in  the  Gulf  of  Mexico.  Characteristics  of  the  near-inertial  currents 
observed  at  the  EFGB  are  expected  to  be  similar  at  the  other  banks. 
In  combination  with  the  mean  currents  and  passing  eddies,  in  the 
vicinity  of  the  banks,  NlCs  probably  play  a  significant  role  in  the 
transport  and  mixing  of  physical  and  biochemical  properties 
around  the  numerous  banks  at  levels  significantly  greater  than 
over  most  of  the  continental  shelf.  Hence,  they  may  have  a 
pronounced  effect  on  the  ecosystem  of  the  Gulf  of  Mexico  shelf 
region.  In  particular,  NlCs  are  most  intense  during  spring  and 
summer,  the  period  of  most  active  biological  activity,  such  as  the 
mass  coral  spawning  event  that  happens  annually  7-10  days  after 
the  full  moon,  typically  in  the  month  of  August. 
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